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Abstract.  It  is  argued  that  the  relativistic  definitions  of  parallax,  proper 
motion  and  radial  velocity  consistent  with  an  accuracy  of  1  fi as  should  be 
considered  only  within  a  well-defined  algorithm  of  relativistic  reduction 
of  observational  data.  Such  an  algorithm  is  formulated  and  the  corre¬ 
sponding  definitions  of  astrometric  parameters  are  discussed. 


1.  Introduction 

The  accuracy  of  future  space-based  astrometric  observations  is  expected  to  at¬ 
tain  a  level  of  1  fi  as.  In  this  paper  we  briefly  describe  a  relativistic  model  of 
space-based  optical  positional  observations  valid  at  such  a  high-accuracy  level 
and  suggest,  in  particular,  definitions  of  parallax,  proper  motion  and  radial 
velocity  compatible  with  general  relativity  at  a  level  of  1  fi  as.  Although  the  def¬ 
initions  are  quite  simple  (see  Klioner  and  Kopeikin,  1992),  their  interpretation 
at  such  a  high  level  of  accuracy  is  rather  tricky.  Parallax  and  proper  motion 
are  no  longer  two  independent  effects.  Second-order  effects  due  to  parallax  and 
proper  motion  as  well  as  the  effects  resulting  from  interaction  between  the  two 
effects  are  important.  Moreover,  parallax,  proper  motion  and  other  astromet¬ 
ric  parameters  are  defined  operationally  and  have  some  meaning  only  within 
the  particular  model  of  relativistic  reductions  chosen.  That  is  why  the  whole 
relativistic  model  of  observations  must  be  considered.  It  is  also  clear  that  in 
order  to  convert  observed  proper  motion  and  radial  velocity  into  true  tangential 
and  radial  velocities  of  the  observed  object,  additional  information  is  required. 
Since  that  information  is  not  always  available,  the  concepts  of  “apparent  proper 
motion”,  “apparent  tangential  velocity”  and  “apparent  radial  velocity”  are  sug¬ 
gested.  These  concepts  represent  useful  information  about  the  observed  object 
and  should  be  distinguished  from  “true  tangential  velocity”  and  “true  radial 
velocity”.  Definitions  of  all  these  concepts  are  discussed  below. 

Throughout  this  paper  the  following  constraints  on  the  various  parameters 
are  used  to  decide  if  a  particular  effect  should  be  retained  to  attain  the  accuracy 
of  1  n as:  (1)  barycentric  position  of  the  observer:  |xs|  <  1  a.u.  +  2  •  106  km;  (2) 
barycentric  velocity  of  the  observer:  |xs|  <  4  •  104  m/s;  (3)  parallax  7r  <  1". 


2.  Relativistic  model  of  positional  observations 

A  relativistic  model  for  optical  positional  observations  having  an  accuracy  of  1 
H as  is  outlined  below.  It  is  assumed  that  the  observations  are  performed  from 
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an  Earth  satellite  or  a  space  station  whose  position  xs  relative  to  the  Barycen- 
tric  Reference  System  (BRS)  of  the  solar  system  is  known  for  any  moment  of 
bary centric  coordinate  time  t  =  TCB.  The  relativistic  model  to  be  discussed 
below  consists  of  several  subsequent  steps  accounting  for  the  following  effects: 

(1)  aberration  (effects  vanishing  together  with  the  barycentric  velocity  of 
the  observer).  This  step  converts  the  observed  direction  to  the  source  s  into  the 
unit  BRS  coordinate  velocity  of  the  light  ray  n  at  the  point  of  observation  xs; 

(2)  gravitational  light  deflection  for  the  source  at  infinity.  This  step  converts 
n  into  the  unit  direction  of  propagation  a  of  the  light  ray  infinitely  far  from  the 
solar  system; 

(3)  coupling  of  finite  distance  to  the  source  and  the  gravitational  light  de¬ 
flection  in  the  gravitational  held  of  the  solar  system.  This  step  converts  a  into 
a  unit  coordinate  BRS  direction  k  going  from  the  source  to  the  observer; 

(4)  parallax.  This  step  converts  k  into  a  unit  BRS  direction  1  going  from 
the  barycenter  of  the  solar  system  to  the  source; 

(5)  proper  motion.  This  step  provides  a  reasonable  parameterization  of  the 
time  dependence  of  1  caused  by  the  motion  of  the  source  relative  to  the  BRS. 

2.1.  Aberration 

The  first  step  is  to  get  rid  of  the  aberrational  effects  related  to  the  BRS  velocity 
of  the  observer.  Let  s  denote  the  unit  direction  toward  the  source  as  observed  by 
the  observer.  Let  p  be  the  BRS  coordinate  velocity  of  the  photon  in  the  point 
of  observation,  (p  is  directed  from  the  source  to  the  observer.)  The  unit  BRS 
coordinate  direction  of  the  light  ray  n  =  p/|p|  can  be  computed  as 

S  =  _n+  “fl  +  ^f2  +  ^f3  +  0{C~A).  (1) 

This  formula  contains  relativistic  aberrational  effects  f \  =  f)(n,xs)  up  to  third- 
order  in  1/c.  Explicit  form  of  f \  can  be  found  in  (Klioner,  1991;  Klioner, 
Kopeikin,  1992).  Because  of  the  first-order  aberrational  terms  (classical  aberra¬ 
tion)  the  BRS  coordinate  velocity  of  the  satellite  must  be  known  to  an  accuracy 
of  10-3  m/s  in  order  to  attain  an  accuracy  of  1  //as.  For  a  satellite  with  the 
BRS  velocity  |xs|  <  40  km/s,  the  first-order  aberration  is  of  the  order  of  28", 
the  second-order  effect  may  amount  to  3.6  mas,  and  the  third-order  effects  are 
~  1  lias.  Note  also  that  the  higher-order  aberrational  effects  are  nonlinear 
with  respect  to  the  velocity  of  the  satellite  and  cannot  be  divided  into  “annual” 
and  “diurnal”  aberrations  as  can  be  done  with  the  first-order  aberration  for  an 
Earth-bound  observer. 

2.2.  Gravitational  light  deflection 

The  next  step  is  to  account  for  the  general-relativistic  gravitational  light  defec¬ 
tion,  that  is  to  convert  n  into  the  corresponding  unit  direction  a  of  the  light 
propagation  infinitely  far  from  the  gravitating  sources.  Relativistic  equations  of 
light  propagation  can  be  written  in  the  form 

xp(Z)  =  Xp(Z0)  +  ca(t  -t0)  +  Sxp(t),  (2) 

where  to  is  the  moment  of  observation,  xp(Zo)  is  the  position  of  the  photon  at  the 
moment  of  observation  (This  position  obviously  coincides  with  the  position  of 
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the  satellite  at  that  moment  xp(/o)  =  xs(Zo).),  is  the  unit  coordinate  direction 
of  the  light  propagation  at  past  null  infinity 

cr  =  lim  c-1  xp(Z),  (3) 

t — *  —  OO 

and  Sxp  is  the  sum  of  all  the  gravitational  effects  in  the  light  propagation 
(Sxp(to)  =  0,  lim  Sxp(t)  =  0).  If  time  dependence  of  the  gravitational  fields 

t — *  —  OO 

produced  by  the  matter  outside  the  solar  system  is  negligible,  those  fields  can  be 
neglected  while  computing  Sxp.  Otherwise  the  external  fields  must  be  taken  into 
account  (e.g.,  for  the  observations  of  close  binary  stars,  where  the  gravitational 
held  of  the  companion  can  cause  an  additional  time-dependent  light  deflection). 

The  largest  contribution  to  Sxp  due  to  the  solar  system  gravity  comes  from 
the  spherically  symmetric  components  of  the  gravitational  fields  of  the  massive 
bodies  (see,  e.g.,  Klioner  and  Kopeikin,  1992).  Depending  on  the  problem  and 
the  required  accuracy,  one  can  additionally  consider  the  effects  of  quadrupole 
moments  of  the  bodies,  their  rotational  and  translational  motion  as  well  as 
the  post-post-Newtonian  effect  (Klioner,  1991a;  Klioner,  1991b;  Klioner  and 
Kopeikin,  1992).  It  is  easy  to  see  that  not  only  the  major  planets  should  be 
taken  into  account  here,  but  also  the  Moon,  Ceres  and  a  dozen  large  satellites  of 
giant  planets.  Detailed  estimates  of  various  effects  will  be  published  elsewhere. 

Coordinate  velocity  of  the  photon  can  be  obtained  by  differentiating  Eq. 
(2):  xp  =  ccr  +  Sxp  and  then  normalized  to  give  the  unit  coordinate  direction  of 
the  light  propagation  at  the  moment  of  observation 

n  =  cr  +  g(n,<Sxp(/0)).  (4) 

2.3.  Coupling  of  finite  distance  and  gravitational  deflection 

The  next  step  is  to  convert  a  into  a  BRS  direction  from  the  source  to  the  ob¬ 
server.  Let  xs(Z)  be  the  coordinate  of  the  satellite  at  the  moment  of  observation 
t  and  X(T)  be  the  position  of  the  source  at  the  moment  of  emission  T  =  T(t) 
of  the  observed  signal.  Let  us  denote 

k(Z)  =  R(Z)/iE(Z),  R(Z)  =  xs(Z)  —  X(T),  R(t)  =  |R(Z)|.  (5) 

Vector  k  is  related  to  cr  as  (Klioner,  1991a): 

a  =  k  —  k  x  (Sxp  x  k)  /  R  +  0(c~4).  (6) 

The  only  effect  in  Sxp  to  be  accounted  for  here  is  the  post-Newtonian  gravita¬ 
tional  deflection  from  the  spherically  symmetric  part  of  the  gravitational  held  of 
the  Sun.  The  explicit  formulas  will  be  published  elsewhere.  The  effect  amounts 
to  10  g, as  for  a  source  situated  at  a  distance  of  1  pc  and  observed  at  the  limb  of 
the  Sun.  One  can  check  that  the  effect  is  larger  than  1  //as  if  |X|  <  10  pc  and 
the  source  is  observed  within  2.3°  from  the  center  of  the  Sun.  If  at  least  one  of 
these  conditions  is  violated  one  can  put  k  =  cr. 

2.4.  Parallax 

Now  we  have  to  get  rid  of  the  parallax  (that  is  to  transform  k  into  a  unit  vector 
1  directed  from  the  bary center  of  the  solar  system  to  the  source): 

1(Z)  =  X(T)/|X(T)|. 


(7) 
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Let  the  parallax  of  the  source  be  defined  as 

7r(t)  =  la.u./|X(T)|,  (8) 

the  parallactic  parameter  II  be  given  by 

n{t)  =  7r(/)  xs(Z)  /  1  a.u.,  (9) 

and  hnally  the  observed  parallactic  shift  of  the  source  be  defined  as 

7 r(Z)  =  1  (/)  X  ( n(t )  X  1(f)).  (10) 

With  these  definitions  to  sufficient  accuracy  one  has 

k  =  -1(1  -  |tt|2  /2)  +  7T  (1  +  1  •  n)  +  0(ir3).  (11) 

The  second-order  effects  in  (11)  proportional  to  7r2  are  less  than  3  fias  if  |X|  >  1 
pc.  The  second-order  terms  can  be  safely  neglected  if  |X|  >2  pc. 

2.5.  Proper  motion 

The  last  step  of  the  algorithm  is  to  provide  a  reasonable  parametrization  of  the 
time  dependence  of  1  and  7r  caused  by  the  motion  of  the  source  relative  to  the 
solar  system  barycenter.  The  following  simple  model  for  the  coordinates  of  the 
source  is  adopted  here 


X(T)  =  Xo  +  VAT  +  (!)(AT2),  (12) 

where  AT  =  T  —  To,  Xo  =  X(To),  and  V  is  the  BRS  velocity  of  the  source 
evaluated  at  the  initial  epoch  To.  This  model  allows  one  to  consider  single 
stars  or  components  of  gravitationally  bounded  systems,  periods  of  which  are 
much  larger  than  the  time  span  of  observations.  In  more  complicated  cases 
special  solutions  for  binary  stars,  etc.  should  be  considered.  Depending  on  the 
source  and  time  span  of  the  observations  higher-order  terms  in  (12)  can  also  be 
considered.  Substituting  (12)  into  definitions  of  1  and  7r  one  gets 

v r(Z)  =  7T0(1  — lo'V/X0AT  +  e>(AT2)),  (13) 

1  =  l0  +  lox(Vxlo)/X0AT  +  e>(AT2),  (14) 

where  Xo  =  |Xo|,  vro  =  7r (ffi)  =  1  a.u./Xo. 

The  signals  emitted  at  moments  To  and  T  are  received  by  the  observer 
at  moments  ffi  and  T  respectively.  Corresponding  moments  of  emission  and 
reception  are  related  by  the  equations 

c(/-T)  =  |xs(/)-X(T)|  (15) 

and  similar  equation  for  ffi  and  To-  The  relativistic  effects  in  (15)  can  be  shown 
to  be  negligible.  Let  us  denote  At  =  t  —  to  as  the  time  span  of  observations 
corresponding  to  the  time  span  of  emission  AT.  These  two  time  intervals  are 
related  as 

At  =  (l  +  c-1 10  •  V)  AT  -  c-1  [xs(Z)  -  xs(Z0)]  •  lo  +  0(AT2).  (16) 
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The  first  term  in  (16)  is  linear  with  respect  to  AT,  the  second  term  represents 
a  quasi-periodic  effect  with  an  amplitude  of  about  500  s,  giving  a  quasi-periodic 
term  in  apparent  proper  motion  of  the  source  (see,  Klioner  and  Kopeikin,  1992). 
Eq.  (16)  results  from  a  Taylor  expansion  of  (15).  Which  terms  of  such  an 
expansion  are  important  depends  on  many  factors.  For  example,  for  a  large 
time  span  of  observations  terms  quadratic  in  AT  may  become  important. 

It  is  easy  to  see  from  (13)  that  time  dependence  of  parallax  7r  can  be  used 
to  determine  radial  velocity  of  the  source.  This  question  has  been  investigated 
in  more  detail  in  Dravins  et  al.  (1999).  The  “true”  tangential  and  radial  com¬ 
ponents  of  barycentric  velocity  V  of  the  source  can  be  defined  by 

vtan  =  lo  x  (V  x  lo),  Erad  =  lo  •  v.  (17) 

Eqs.  (13)-(14)  can  be  combined  with  (16)  to  get  the  time  dependence  of  1  and 
7T  as  seen  by  the  observer.  Collecting  terms  linear  with  respect  to  At  we  get  the 
definition  of  apparent  tangential  velocity  V(^n  as  the  linear  term  in  1(f),  and  of 
apparent  radial  velocity  Wad  as  the  linear  term  in  vr(t): 


VtTn  =  Vtan  (l  +  C"1  Wad)  ,  Wad  =  ^rad  (l  +  C"1  Wrad)  •  (18) 

With  these  definition  the  simplest  models  for  7 r(/)  and  1(f)  as  seen  by  the  observer 
read  (the  higher-order  terms  are  neglected  here): 

7T(t)  =  vr0  -  T^l  (Wad  /  1  a.e.)  At  +  . . . ,  (19) 

1(1)  =  lo  +  Map  At  +  MapC-1  (tXs(t)  -  XS(to)]  '  lo)  +  •  •  •,  (20) 

Map  =  ^0  (VtaW/l  a.e.).  (21) 

Apparent  proper  motion  is  denoted  fi  here.  The  factor  (l  +  c-1  Wad)  1  in  (18) 


has  been  discussed  in,  e.g.,  Stumpff  (1985)  and  Klioner  and  Kopeikin  (1992). 
This  factor  may  become  very  large  and  is  one  of  the  possible  explanations  of 
the  well-known  problem  of  superluminal  motions  in  quasars  and  active  nuclei  of 
galaxies. 

If  both  and  Wad  can  be  determined  from  observations  one  can  imme¬ 
diately  restore  the  “true”  velocities  Vtan  and  Wad-  However,  even  if  it  is  not 
the  case  and  Wad  are  useful  by  themselves.  Note  that  the  radial  velocities 
as  measured  by  Doppler  observations  are  affected  by  a  number  of  factors  not 
influencing  positional  observations  and  do  not  coincide  with  either  Vrad  or  Wad- 
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